Hamiltonian of a spinning test-particle in curved spacetime 
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Using a Legendre transformation, we compute the unconstrained Hamiltonian of a spinning test- 
particle in a curved spacetime at linear order in the particle spin. The equations of motion of this 
unconstrained Hamiltonian coincide with the Mathisson-Papapetrou-Pirani equations. We then use 
the formalism of Dirac brackets to derive the constrained Hamiltonian and the corresponding phase- 
space algebra in the Newton- Wigner spin supplementary condition (SSC), suitably generalized to 
curved spacetime, and find that the phase-space algebra (q, p, S) is canonical at linear order in 
the particle spin. We provide explicit expressions for this Hamiltonian in a spherically symmetric 
spacetime, both in isotropic and spherical coordinates, and in the Kerr spacetime in Boyer-Lindquist 
coordinates. Furthermore, we find that our Hamiltonian, when expanded in Post-Newtonian (PN) 
orders, agrees with the Arnowitt-Deser-Misner (ADM) canonical Hamiltonian computed in PN 
theory in the test-particle limit. Notably, we recover the known spin-orbit couplings through 2.5PN 
order and the spin-spin couplings of type Skch S (and S^err) through 3PN order, 5*Kerr being the 
spin of the Kerr spacetime. Our method allows one to compute the PN Hamiltonian at any order, 
in the test-particle limit and at linear order in the particle spin. As an application we compute it 
at 3.5PN order. 

PACS numbers: 04.25.D-, 04.25.dg, 04.25.Nx, 04.30.-w 
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I. INTRODUCTION 

The dynamics of spinning bodies in general relativity 
is a complicated problem which has been investigated 
in several papers during the last seventy years, start- 
ing from the pioneering work by Mathisson [l|, Papa- 
petrou [IS II, Pirani [|, Tulczyjew @,0| and Dixon]!]. 
Spin effects on the free motion of a test particle were first 
derived in the form of a coupling to the spacetime curva- 
ture in Refs. @, [!,[!]• The computation assumes that the 
test-particle can be described by a polc-dipolc energy- 
momentum tensor [!, @] , thus neglecting the quadrupole 
moment (and higher multipole moments) and providing 
spin couplings only at linear order in the test-particle's 
spin. 

The two-body dynamics of spinning objects can also 
be computed in post-Newtonian (PN) theory [9[, which 
is basically an expansion in powers of v/c and GM / (c 2 r), 
where v is the characteristic velocity of the system and 
r is the binary's separation. Currently, spin couplings 
have been computed in the two-body equations of mo- 
tion through 2.5PN order 0, EL El 03, 0, d, and 
in the Arnowitt-Deser-Misner (ADM) canonical Hamilto- 
nian through 3PN order [11 [U, Gl, HO] and partially 
at higher PN orders [2l|, |22| |. These coupling terms agree 
with those computed via effective- field-theory techniques 
at 1.5PN, 2PN and 3PN order [H, [H [H, Hi ■ 

The main motivation for describing as accurately as 
possible the dynamics of a binary system of spinning com- 
pact bodies in general relativity comes from the forth- 
coming observation of gravitational waves with ground 
and space-based detectors. In particular, LIGO, Virgo 
and GEO could observe signals emitted by stellar-mass 
black-hole and neutron-star binaries, and LISA could de- 
tect signals from supcrmassivc black-hole binaries and 



extreme-mass ratio binaries. 

In this paper we compute the Hamiltonian of a test- 
particle in a curved background spacetime, including all 
couplings linear in the test-particle's spin. Starting from 
the Lagrangian given in Ref. (27[, we apply a Legendre 
transformation to derive the unconstrained Hamiltonian. 
The Hamiltonian is unconstrained in the sense that the 
test-particle's spin variables are given by an antisym- 
metric tensor S^, which a priori contains six degrees 
of freedom instead of three. It is well-known that in or- 
der to fix the unphysical degrees of freedom associated 
with the arbitrariness in the definition of , a choice 
must be made for the so-called spin supplementary con- 
dition (SSC). The arbitrariness can be interpreted, in the 
case of extended bodies 1 , as the freedom of choosing the 
point, internal to the body, whose motion is followed [29| . 

Building on the work by Hanson and Regge [3(| and 
generalizing the Newton- Wigner (NW) SSC to curved 
spacetime, we then derive the constrained Hamiltonian 
and the corresponding Dirac brackets, which should re- 
place the Poisson brackets when computing the equations 
of motion from that Hamiltonian. Quite interestingly, we 
find that the NW SSC leads, at least at linear order in 
the particle spin, to canonical Dirac brackets, i.e. the 
standard sympletic structure for a set of dynamical vari- 
ables (q,p, S). As a consistency check of our results we 
also compare our constrained Hamiltonian with the ADM 



1 It should be stressed that any spinning "particle" must actually 
have a small non-finite size. An intuitive argument for this can 
be found in Ref. [2H , Ex. 5.6, where it is shown that any spinning 
body must have a minimal size in order not to rotate at velocities 
larger than c. A more rigorous proof can be found in Ref. l2ll , 
Sec. 2. 



2 



canonical Hamiltonian for spinning bodies, as computed 
in PN theory through 3PN order. In addition we provide 
explicit expressions for the Hamiltonian of a spinning 
particle moving in a generic spherically symmetric space- 
time (using both isotropic and spherical coordinates), as 
well as in the Kerr spacetime (in Boyer-Lindquist coor- 
dinates). 

Another important application of this work will be de- 
veloped in a subsequent paper, where we will use the 
Hamiltonian derived here to build a new effective-one- 
body Hamiltonian [3l|, [H, [H, HH for spinning objects. 
This application is crucial to take full advantage of the 
analytical and numerical treatment of the dynamics of 
spinning bodies throughout the inspiral, merger and ring- 
down, and build accurate templates for the search of 
gravitational waves with ground-based and space-based 
detectors. 

The paper is organized as follows. In Sec. [Til we 
briefly summarize our notations. In Sec. IIIII we ap- 
ply a Legendrc transformation to compute the uncon- 
strained Hamiltonian and show that the equations of mo- 
tion that follow from it coincide with the well-known 
Mathisson-Papapetrou-Pirani (MPP) equations of mo- 
tion. In Sec. IIV1 after reviewing the Dirac bracket 
formalism, we derive the constrained Hamiltonian and 
the corresponding Dirac brackets using the generalized 
NW SSC. In Sec. |V] we specialize our results to spher- 
ically symmetric spacetimes and to the Kerr spacetime 
in Boyer-Lindquist coordinates. In Sec. IVII we restrict 
ourselves to the Kerr spacetime in ADM coordinates, ex- 
pand the Hamiltonian computed in the NW SSC in a 
PN series through 3.5PN order and find agreement with 
the ADM canonical Hamiltonian in the test-particle limit 
through 3PN order. Section IVIII summarizes our main 
conclusions. 
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e- A e v B g^ v = i]ab , (2.1) 

where t\tt = —1, Vti = 0, J]u = 5u (Su being the 
Kroncckcr symbol). Thus the internal tetrad space is 
Lorentz invariant, i.e. one can obtain any tetrad from an 
existing one by applying a Lorentz transformation e' A = 
A A B eB, where 

A A C A C B =A C A A B C = 8 B . (2.2) 

Internal tetrad indices denoted with the uppercase Latin 
letters A, B, C and D always run from to 3, while 
internal tetrad indices with the uppercase Latin letters /, 
J, K and L, associated with the spacelike tetrad vectors, 
run from 1 to 3 only. The timelike tetrad index is denoted 
by T. 

Tetrad indices are raised and lowered with the metric 
Vab [e.g., e A = t\ab (e s ) M ]. With this convention the 
relation (|2.ip can be easily shown to be equivalent to the 
completeness relation 

= K ■ (2.3) 

We will denote the projections of a vector V onto the 
tetrad with V A = and similarly for tensors of 

higher rank, as well as Christoffcl symbols. Partial 
derivatives will be denoted with a comma or with <9, co- 
variant derivatives with a semicolon, while total covariant 
derivatives with respect to a parameter a will be denoted 
by D/Da. Finally, we will denote the operation of an- 
tisymmctrization with respect to the indices jj, and v as 
A " 11- = :.l " />"' -A-i*B u ~)/2. 

We use geometric units G = c = 1 throughout the pa- 
per, except in Sec. I VII where the factors of c are restored, 
playing the role of PN book-keeping parameters. 

III. UNCONSTRAINED HAMILTONIAN 



In this section we derive the unconstrained Hamilto- 
nian by applying a Legendre transformation to the La- 
grangian describing the motion of a spinning particle in 
a generic curved spacetime. 



II. NOTATIONS 



Throughout this paper, we will use the signature 
(—,+,+,+) for the metric. Spacetime tensor indices 
(ranging from to 3) will be denoted with Greek let- 
ters, while spatial tensor indices (ranging from 1 to 3) 
will be denoted with lowercase Latin letters. Also, we 
will often use t as alternate for the timelike index 0. 

We define a tetrad field as a set consisting of a timclike 
future-oriented vector ij, and three spacelike vectors ij 
(I = 1, 3) — collectively denoted as e A (A = 0, ...,3) 



A. The Lagrangian and the 
Mathisson-Papapetrou-Pirani equations 

Building on the classic work of Hanson and Regge [3(3| 
which analyzes the dynamics of a relativistic top in a flat 



2 We use the notation e A to denote any choice of tetrad given a 
background spacetime. The tetrad without the tilde e A refers 
to a special tetrad, namely the one carried by the test particle. 
The tetrad e A is special in the sense that it is a dynamical vari- 
able whose evolution along the worldline is prescribed by some 
Lagrangian. 
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spacctime, Porto showed in Ref. [27J that the equations 
of motion of a spinning particle in curved spacetime can 
be obtained from the action 



S = J L(ai, a 2 , as, a 4 ) da , 



(3.1) 



a being a parameter along the representative worldlinc. 
The Langrangian L is a function of the four Lorentz- 
invariant scalars 



«i 

2 
0.3 
0,4 



(3.2) 
(3.3) 
(3.4) 
(3.5) 



where = dx^/da is the tangent vector to the repre- 
sentative worldline, and where the antisymmetric tensor 



W describes how the tetrad 
rotates along the worldline: 



carried by the particle 



AB u De B 

V €a Da 



da 



a/3 9 



(3.6) 



Moreover, the action (|3.1| is assumed to be 
reparametrization-invariant (i.e. its form must be 
independent of the particular parameter used to follow 
the particle's worldline), which translates in the require- 
ment that the Lagrangian L be a homogeneous function 
of degree one in the "velocities" and Vt^ [3(| ■ Porto 
then shows that if one defines the four-momentum 
vector and the spin tensor of the particle as 3 



dL 



Sau — 2 



dL 



(3.7) 
(3.8) 



[note that p^ is not the momentum conjugate to the co- 
ordinates because f2 M " depends on u M , as can be seen 
in Eq. (|3.6|) ]. then a variation of the action with respect 
to which preserves the defining property (|2.1[) of a 
tetrad gives the precession equation for the spin tensor 



Da 



5" 



p^ u — p u H 



(3.9) 



Thus one precisely recovers the well-known MPP equa- 
tions from the action (|3.ip . which therefore encodes the 
dynamics of a spinning test-particle in curved spacetime, 
at linear order in the particle's spin. 

Notice however that the set of Eqs. (|3.9|) - (|3.10p con- 
sists of ten equations and thirteen independent variables 
(p^, and , subject to the normalization constraint 4 
of the tangent vector u M ) and is therefore not closed. This 
undcrdetermination can be addressed by imposing a SSC, 
which is typically expressed as 



S^LOy = 0, 



(3.11) 



where uj v is some suitably chosen timelike vector. Equa- 
tion (|3.11|) contains three independent constraints, and 
is therefore expected to reduce the number of indepen- 
dent variables from 13 to 10, thus closing the system of 
Eqs. (|3.9[) - (|3.10p . This is indeed what happens, as the 
requirement that Eq. (13.111) be valid at all points along 
the worldline implies the following implicit relationship 
between p^ and 



p» = 



1 



OJyU" 



{uj vV v )u» - 



Du 



Da 



(3.12) 



It should be stressed once again that it is the underde- 
termination of the unconstrained MPP system that al- 
lows one to impose any constraint of the form (|3.1ip . 
and that the constraint will be automatically conserved 
by the time evolution of the system because of Eq. (|3.12l) . 
Of course different constraints of the form (|3.1ip will pro- 
duce different systems of equations describing the evolu- 
tion of the particle's worldlinc. The physical reason for 
this is easy to understand: the SSC (|3.1ip binds the test- 
particle described by the Lagrangian to a specific, SSC- 
dependent, worldline lying inside the worldtube spanned 
by the spinning body, namely the center of energy of the 
body as seen by an observer with four- velocity parallel to 
lo^ (see e.g. Ref. [29[ for a lucid discussion of the physical 
meaning of SSCs). 



B. Deriving the Hamiltonian through a Legendre 
transformation 



The second equality in Eq. (|3.9p follows from defini- 
tions (|3.7p and (|3.8p . and from the fact that the La- 
grangian in the action (|3.1|) depends only on a\, a 2 , 03 
and 04 (30j |. Moreover, a variation of the action with 
respect to the particle's position x^ gives [2?| 



Dpt* 
~Da~ 



l -R\ M u«S^ 



(3.10) 



It is convenient to rewrite the action (|3.1|) as 



S = I L[ x^,^,-^- ) da, 



(3.13) 



where the Langrangian L can be now considered as a 
function of the coordinates the four- vector u M = 



Because of reparametrization invariance of the action (|3.1| l. these 
definitions maintain the same form whatever parameter a is cho- 
sen along the wordline, as appropriate for physical quantities like 
the four-momentum and the spin. 



4 For example one is free to select a parameter a = r such that 
u^u^ = —1, since the action 113.11 1 is reparametrization invariant. 
Any other choice of parameter simply yields a different normal- 
ization constraint u^u^ = — (da/dr) -2 . 
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dx^/da, the six parameters <j) a and their time deriva- 
tives. The set {4> a } consists simply of the parameters of 
the internal Lorentz transformation describing the orien- 
tation of the tetrad field e A carried by the particle with 
respect to an arbitrary, but fixed, reference tetrad field 
e A (x) covering the whole spacetime 5 . Therefore, the 
tetrad carried by the particle is given by 



[(fax) 



A A B (fa£ B (x), 



(3.14) 



where A A B is a Lorentz transformation. We also note 
that the parameters <p a and their time derivatives enter 
the Lagrangian only through the antisymmetric tensor 
fJA" / ] which we write explicitly as 



fi<- = V AB e>X(fax) 



b a de v B 



^(fax)+u e B>0 (fax) 



da 



(3.15) 



To construct the Hamiltonian we need to choose a par- 
ticular 3+1 decomposition of the background metric. We 
take a = t, where t is the time coordinate of that partic- 
ular decomposition. Using reparametrization invariance, 
we can write 



S = I L [ a*,*,", ) da, 



L(x\v\cj) a ^ a ,t)dt, 



(3.16) 



where . 



= t, u° = 1, u l = v l = dx l /dt and <j) a = dfa-jdt. 
The configuration space of the spinning particle therefore 
consists of the set {x l , cj) a }. The total variation of the 
Lagrangian considered as function of x l , v l , <fi a and (j> a is 



. T dL ■ dL ■ dL , dL ■ „ 

bt Fit 
= Q-i ^ + p i 5vi + gTa *P + . (3-17) 



where we denoted by Pi and P^ the momenta conjugate 
to x % and cf> a , respectively. The total variation of the 
Lagrangian considered as function of x i , v l and W is 
instead 



Using Eq. (|3.15[) . Eq. (|3.18|) can be rewritten as 



SL 



dL 
dx l 
dL_ 
dv l 
dL 



dL 



dL 



dW v 



dL 



n BQ^ V 
dW v 



x,v dx % 



,v dv 1 



a;, 0,0 



5x l 



x,v d(j) a 
dW v 



X,V,(f) 



(3.19) 



Comparing Eq. (f3TT7|) with Eq. (|3~19)) . and using 
Eqs. (|3.7[) . (|3.8p and Eq. (|3.15|) . we obtain the conjugate 
momenta 



P 



1 

- n AB <? p v 
2 '/ Bpv <-A C B;i I 



(3.20) 



and 



2 V -V e A 



1 



Q \AB ~H ~v 
2 A a L A C B i 

where we have introduced the tensor 

1 ~yl ~B 

E\^ v = - t]ab e p e v . x , 



(3.21) 



(3.22) 



which is antisymmetric in the last two indices, and the 
antisymmetric tensor [30l ] 



g A CB 



(3.23) 



A necessary condition to go from the Lagrangian for- 
malism to the Hamiltonian one in the usual way (i.e. by 
means of a Legendre transformation) is that the Lan- 
grangian is regular [35| . i.e. it satisfies 6 



det ( ) 



(3.24) 



where q = (x l , 4> a ). Under this condition, we can perform 
the usual Legendre transformation to get the Hamilto- 
nian 



H = P, v l + P 6 



L. 



(3.25) 



SL 



dL 

dx l 



n dv 1 



5v l 



dL 



dQ^ 



(3.18) 



5 In what follows we will prove that the equations of motion are 
independent of the choice of this tetrad field. This had to be 
expected, based on Refs. [27l, l3dH . 



° While this condition is sufficiently generic to leave our La- 
grangian essentially undetermined, it should be noticed that 
there are famous examples in physics where this regularity con- 
dition does not hold, such as the electromagnetic field (see for 
instance Ref. [36l |. chapter 5), the Dirac field (see for instance 
Ref. [37J, problem 9. 2d), the Schrodinger equation (see Ref. [3H 
and references therein) and general relativity (see for instance 
Ref. [H, chapter 9). 
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Since L is homogeneous of degree one in the "velocities" 
[because of the reparametrization invariance of the ac- 
tion (|3.1|) ]. Euler's theorem implies that 



dL 
duf* 



dL 



(3.26) 



where we have used the definitions (|3.7p and (|3.8p , as 
well as the fact that with our time-slicing u = 1 and 
u l = v\ Using now Eqs. (j3~2"6l) . (|3"^Uj) (|3~2Tj ) and (|3~T5l) 
(with u° = 1 and u l = z/) in Eq. (|3.25p . simple algebra 
allows one to write the Hamiltonian as 



tt 1 AB o a 3 

H = -p t - -v S a(3 e A e^, , 



-Pt- ^l AB S al3 ~e a A e 
-p t - Ek^S" 



a -a 

B;t ' 



(3.27a) 

(3.27b) 
(3.27c) 



where the covariant derivative with respect to t in the 
second term of Eq. (|3.27a[) above is a shorthand for co- 
variant derivative with respect to x° = t, i.e. one can 
pull the Lorentz transformation A A B ((f>) outside the co- 
variant derivative as it is independent of x°. It should be 
noted that using the tensor Ex^ defined in Eq. (|3 . 22[) . 
one can combine H and Pj into a four-vector P Q such 
that 



P a = {-H, Pi) = Pa 



TP QV-V 



(3.28) 



The MPP equations of motion can be derived from the 
Hamiltonian (|3.27c|) as follows. On one hand we have 



dry 

dt 



dH dL 



dL 



X,U,tp 



dW 



(3.29) 



where the second equality follows from the definition 
of the Hamiltonian (|3.25[) with the regularity condi- 
tion (|3.24p . (One could also derive the second equality by 
comparing the Hamiltonian and Lagrange equations, but 
it should be stressed that these two sets of equations are 
equivalent only if the regularity condition (|3.24|) is satis- 
fied dl|.) Using then Eqs. (|3"fl5j) and ([3~2"T]) . as well as 
the definition (|3.8p . a straightforward computation gives 
the precession equation 



DS^ 
Dt 



(3.30) 



The translational equations of motion can be obtained 
following a similar procedure. In the neighborhood of any 
event located on the particle's worldline wc can choose 
Riemann normal coordinates and write 



<m 

dt 



dpi 
dt 
dH 

dx l 
dL 



dL 

dx l 



dx l 



(3.31) 



where the last equality follows from the compatibility of 
the metric with the connection, i.e. g^-i = 0, which be- 
comes g^u.i = in Riemann normal coordinates. 7 Mak- 
ing use of Eq. (|3.15p and using the fact that in Riemann 
normal coordinates = 0, while their derivatives are 
non-zero, we get 



dp t 
dt 



2 Pict/37 u ^ ^ , 



(3.32) 



where the Riemann tensor term arises from the deriva- 
tives of the Christoffcl symbols appearing in Eq. (13.151) . 
Rewriting Eq. (|3 . 32[) in a generic coordinates system, we 
immediately get the spatial part of the translational MPP 
equations 



Dpi 
Dt 



Q®1 



(3.33) 



The unconstrained equation of motion for p t is obtained 
as follows. One starts from the formal expression 



— = {PU H} + — . 



(3.34) 



In Riemann normal coordinates, the left-hand side is 
equal to Dp t /Dt. To evaluate the right-hand side, 
one makes use of Eq. (|3.27cp to eliminate pt in fa- 
vor of the Hamiltonian and other quantities whose ex- 
plicit expressions in terms of the phase-space variables 
{x\ Pi, 4> a , P,),*} arc known. Straightforward algebra 
then yields 



(3.35) 



which can be combined with Eqs. (|3 . 33[) in the well- 
known equation translational MPP equations 



(3.36) 



Before concluding this section, we provide explicit expres- 
sions for the Poisson brackets of the variables x l ,Pi, 
and A^" = k AB e A e u B . Using the definition of Poisson 
bracket, 



{f,g}^ 



df dg dg df 
3q d-rr <9q div 



(3.37) 



where q = (x 1 , <fi a ) and 7r = (Pj, P^), we trivially have 

{x\Pj} = 5), (3.38a) 
{x\x j } = {P,P,} = 0. (3.38b) 



7 We stress that one is allowed to set g^v^i = in this equation 
as we do not need to take derivatives of it (in which case, of 
course, the terms containing g^ Ut i would give a contribution, as 
in general g^v^j ^ even in Riemann normal coordinates). 
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To compute the Poisson brackets involving let us 
first invert Eq. ([3^21]) [Ij: 



nap _ ~a ~f3 AB p 
a — e A e B Pa r 4> 



where p AB ((j>) satisfies 



Pa A bAB — 2<5 a fc , 
x AB CD _ AC BD AD BC 
A a Pa — V V ~ V V 

Using these relations together with the identity 



(3.39) 

(3.40) 
(3.41) 



UA a UA b _ \AC \ B \AC \ B in 
~Jh& Thha~ ~ A « Af, C ~ A >> Q C- > \ 6Al ) 



which can be immediately derived [301 ] by taking the 
derivative of Eq. (|3.23p . it is straightforward to prove 
that p AB is a realization of the Lie algebra of the Lorentz 
group: 

dp C a° CD ®Pa _ „ACBD RD AC 
Pb ~FUtr~ Pb ~ 'Pa V -Pa V 



Simple algebra then yields 



+p£ d v bc + p* c v ad . 



(3.43) 



{S^,S a0 } = S" a g v& + S v(} g" a 

_ s „p g ua _ S u a ^ (3.44a) 



while using Eqs. (|2.3[) and (|3.39[) we easily obtain 

{S^,Pi} = S a »e£e>X. + S>"*e£eX i , (3.44b) 



{S^,x 1 } = 0. 



(3.44c) 



Finally, it is straightforward to show that A AB satis- 
fies m 

{A AB , x 1 } = {A AB , Pi} = {A AB , A CD } = , (3.45) 
{A AB ,S CD } = A AC V BD -A AD V BC , (3.46) 

or, in terms of A**" = A AB e\e v B 

{A"", x 1 } = {A» u , A af3 } = , (3.47) 
{A"", = A a » ~e A + A» a e A ~e\ t , (3.48) 
{A"", S af3 } = A^g" 13 ~ A"V° ■ (3-49) 

IV. CONSTRAINED HAMILTONIAN 

A. Imposing constraints in phase-space: a 
Dirac-bracket primer 

Let us briefly recall how constraints are imposed in 
the Hamiltonian formalism (a very detailed review on 
the subject can be found in Ref. [391]). Let us consider a 
Hamiltonian -ff(q, 7r,t) living in a 2n-dimcnsional phase 



space and a binary "bracket" operation {..., ...} which is 
antisymmetric, bilinear, and which satisfies the Leibniz 
rule, as well as the Jacobi identity, i.e. 

{A, B} = -{B,A}, (4.1a) 

{aA + bB,C} = a{A,C} + b{B,C}, (4.1b) 

{AB,C} = {A, C}B + A{B, C} , (4.1c) 

and 

{A, {B, C}} + {B, {C, A}} + {C, {A, B}} = . (4.1d) 

In Eqs. (|4.ip . A, B and C are arbitrary phase-space func- 
tions, while a and b are constants. Let us also assume 
that the bracket operation gives the equations of motion 
for a generic phase-space function A through the Hamil- 
ton equations 



dA dA r , m 

If we consider now a set of constraints = 0, i = 1, 
(with m < n) such that the matrix 



(4.2) 



. 2m 



(4.3) 



is not singular 8 , these constraints can be imposed simply 
by replacing the original brackets with the so-called Dime 
brackets. The Dirac brackets are in essence the projec- 
tion of the original symplcctic structure onto the phase- 
space surface defined by the constraints. For two arbi- 
trary phase-space functions A and B, the Dirac brackets 
are given by 



{A, £}db = {A B} + {A, ej {B, 0} [C-% 



(4.4) 



It can be shown (see e.g. Sees. 1.3.2, 1.3.3, and Ex. 1.12 
in Ref. [3!|), that the Dirac brackets are bilinear, an- 
tisymmetric, that they satisfy the Leibniz rule and the 
Jacobi identity, and that they provide the correct equa- 
tions of motion for the constrained system through the 
Hamilton equations 



dA dA XA & _ 



(4.5) 



where A is an arbitrary phase-space function, and where 
the new Hamiltonian H is obtained simply by inserting 
the constraints in the original Hamiltonian H. 

In summary, given a Hamiltonian H and a bracket 
operation (e.g., the Poisson brackets in the case of an 
unconstrained Hamiltonian), in order to impose a set 
of constraints satisfying det(C^) ^ [with C given by 
Eq. (|4.3|) ] . we need to replace the original bracket oper- 
ation with the Dirac bracket operation (|4.4p . and insert 
the constraints directly in the original Hamiltonian. 



8 In the literature, constraints satisfying this condition arc known 
as second class constraints [331 . 
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In Sees. IIVBI and IIV CI we start from the uncon- 
strained Hamiltonian (|3.27cp and the unconstrained al- 
gebra (I338ajl . JUSbJ, (IQ4a|l . (l3~44bl) . (pOIcjl . lpTi7jl . 
(|3.48p . (|3.49p , and use the procedure outlined in this 
subsection to impose the generalized NW SSC. In par- 
ticular, in Sec. IIVBI we compute the Dirac brackets in 
the NW SSC, showing that they are canonical (i.e., they 
reduce to the usual Poisson brackets) at linear order in 
the particle's spin, while in Sec. IIV O we explicitly write 
the constrained Hamiltonian. 



B. Dirac brackets in the generalized 
Newton- Wigner spin supplementary condition 

In this section, we consider the NW SSC generalized 
to curved spacctimc, 

V" = uj v = 0, (4.6) 

with 

t*V=?V-meJ, (4.7) 

where m = \/ —p^p^ is a function of phase space variables 
that we define as the mass of the particle 9 . We stress that 
the vector u) is the sum of two timclikc future-oriented 
vectors and is therefore timelike itself, which implies that 
Eqs. gU) and gj]) do indeed yield a legitimate SSC 
(We recall that with our notation one has e T = —ex, and 
that er is future oriented.) 

While the NW SSC is well-known to be the only SSC 
condition which yields canonical variables in flat space- 
time 10 H3,[13,lll| i there is no a priori guarantee that this 
is the case in curved spacetime. In this section we show 
that the NW SSC does indeed yield canonical variables 
at linear order in the particle's spin. 

Because V* 1 uj^ = 0, only three of the four con- 
straints (|4.6p are independent. Since w is a timelike vec- 
tor, it is natural to take the three independent constraints 
to be the spatial components V % . The constraints V % may 
be viewed as constraints on the momenta P<p°- , as there is 
a one-to-one mapping between the spin tensor S^ v and 
the six momenta conjugate to the 4> a, s. This implies that 
by themselves, the constraints V % do not form a consis- 
tent set of constraints on phase-space: an additional set 
of three constraints must be imposed on the configuration 
coordinates <f> a themselves in order to retain a symplcctic 
structure, i.e. that the constraint hypersurface contains 
the same number of configuration coordinates and con- 
jugate momenta. The additional constraints we choose 



9 Note that at this stage there is no guarantee that this function 
on phase space is a constant of motion. We will show later that 
it is indeed the case, but we emphasize that this is a non-trivial 
result. 

We note that in quantum mechanics and fiat spacetime the NW 
SSC holds a special place [iollilll. 



to impose are given by [23|, [3(| 

X, = (er)„ - - 

= A T A (e A )„-^ = 0. (4.8) 

It is worth pointing out once again that the mass m is a 
function on phase space, and therefore its Poisson brack- 
ets with coordinates and momenta arc non-vanishing. It 
will acquire a special status as a constant of motion (at 
linear order in spin) only at the end of this subsection. 
Equation (|4.8| may be alternatively rewritten as 



A T A = -p M (e A y = ^, (4.9) 



which shows explicitly that it constrains the three ve- 
locity parameters, say </> 4,5 ' 6 , of the Lorentz transforma- 
tion that relates the tetrad carried by the particle to the 
background tetrad. Since A T T is fully determined by 
A T 7 , only three of the four constraints given in Eqs. (|4.8p 
or (|4.9[) are independent 11 . We will take the spatial com- 
ponents = as our three independent constraints on 
the coordinates <\> a . 

In summary, for the generalized NW SSC, the vector 
of constraints is 



i= (V\V 2 7 V 3 , X i,X2,X3)- (4.10) 



In principle, the computation of the matrix C defined 
in Eq. (|4.3p can be performed directly using the un- 
constrained symplectic algebra (|3.38a|) . (|3.38b[) . (|3.44ap . 
([3~44bl) . (j5^4cj) . (|3TT7|) . ([3"^]) . However, since 

the constraints are formulated in terms of the momen- 
tum four-vector p^ rather than the conjugate momenta 
Pi and the Hamiltonian H, it turns out to be quite useful 
to first compute Poisson brackets between p^ and other 
phase space quantities, and then make use of these re- 
sults to compute the matrix C. The relevant Poisson 



1 One can also see this from the fact that x is orthogonal to the 
timclike vector e.T + p/ (mc). Hence only its three spacelike com- 
ponents are independent. 
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brackets are 



{x , Pj } = 






(4.11a) 


{x\Pt} = 


-v\ 




(4.11b) 


{a; 1 , to} = 


_I(p*_ p V 

TO 


• 


(4.11c) 


{Pi, Pj] = 


— 1/?.. <?'"' 




(4.11d) 


{Pi > Pt } = 


2 M 




, (4.11e) 


\Pii 1,1 f — 


-—75 r P ('t) A: - 
TO 


P u ) i 


C4 1 1f"l 


{p t ,m} = 


1 „ i. 




(4.11g) 


{Pi,( e r)j} = 


-r&(er) M , 




(4.11h) 


{Pt,( e T)j} = 


- (-n jvP v + 

771 V 




, (4.Hi) 


{S^, Pi } = 


2S A ^T^ , 




(4.Hj) 


is*",*} = 


-2p [ ^ ] - 25 


A[ Mr ^] k 


(4.11k) 


{s^M = 






(4.111) 




2d] e T J , 




(4.11m) 


{(e T ) l7 TO} = 




+ 










(4.11n) 


{{e T YM = 


--(e T r k ( P k 

TO 




(4.11o) 



The total time derivative dA/dt is then evaluated with 
the help of the unconstrained equations of motion. The 
Poisson brackets (f47TT|) along with Eqs. ([3~1"2|) and (|3.44ap 
yield 



{V*,V j } 



TO 

+0(S 2 ) 
1 / 



(4.13) 



{Xi,Xj} = ^^{piRj\^ -p ] Ri\^)p x S tny 



l 



2m 



m 
(4.14) 



(4.15) 



The remainders scaling as the square of the particle's spin 
are dropped, since the pole-dipole particle model is valid 
only at linear order in the particle's spin. The matrix C 
defined in Eq. (|4.3[) is therefore given by 



C = K + S + 0{S 2 ) 
where the matrices K and S arc defined as 



K = 



3 Q 
Q T 3 



(4.16) 



(4.17) 



where the Poisson bracket between an arbitrary phase 
space function A and the quantity p t is obtained as fol- 
lows 

{A,pt} = {A,-H -lr, AB (e A Ue B ) ., t S^} , 

3A dA 1 AB r A /- \ (~ \ q/3i 



with 



(4.12) and 



(4.18) 



1 

'2m 5 



RkXuvS^ 



(4.19) 



The inverse matrix C 1 can be easily computed at linear 
order in the spin, the result being 

C- 1 = K 1 - K 1 + 0(S 2 ) , (4.20) 

where 

-V 1 )- < 4 - 21 » 



with 

lQ- 1 ] i i = J ^U + ^) ■ (4-22) 

To compute the Dirac brackets between two phase space 
functions, one also needs the Poisson brackets between 
those phase space functions and the constraints. For our 
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purposes, the relevant brackets are given by 



{x\V j } 



-S ij - S 3 %+0(S 2 ), 



5} + ^(p*-p t v?) +<D(S*), 



(4.23a) 



(4.23b) 



{S AB ,V 1 } 
{S AB ^} 



S Al u B ~ S Bi lo a + 0{S 2 ) , (4.23e) 



2 \A~B\ 2Pi [A B\ 

— p l e, H ^ p l u 1 

--s^ A rfUo(s 2 ). 



(4.23f) 



{P t , Vi} = p A S ja e£ ti - S^Tt^ + 0(S 2 ) 



gjk _ gjt 



p 



0{S 2 ) , (4.23c) 



{Pi,Xj} = 



-Lp A ef i +'^p^x 



2m 2 ^ 



1 

m 



(4.23d) 



The matrix (j4~20)) and Eqs. (|4~23| . together with the un- 
constrained algebra given by Eqs. ()3.38j) and (|3.44j) . is 
all one needs to compute the Dirac brackets according to 
Eq. (|4.4p . Our results for the Dirac brackets involving x l 
and Pj are given by 



{x\x j } D B 
{x\Pj}bb 



2p v u v 



{p°u a y 



gij _ git? + gjtP ) + , S 2, = Q ,gtS 
pt pt I 



«5 



s lk - s 



up 



s 



ktp_ 

p 1 



2p v u v 



kl QktP | gltP 



s Kt - s 



(p^aY 
UJ^UJa - 2p V UJ v 



p a ~e a A e A J +0{S 2 )=5)+0{S z ), 



(p a uj a y 



P a e a A ^P^F, 3 + 0(S 2 ) = 0(S 2 ) 



(4.24a) 
(4.24b) 
(4.24c) 



The crucial point now is that Eq. (|4.7p implies lo^oj^ = 
2p M <x> A i, and therefore all terms linear in the particle's spin 
on the right-hand side of Eqs. (|4.24p vanish. Hence the 
Dirac bracket algebra between x 1 and Pj is canonical up 
to terms quadratic in the particle's spin. 



The Dirac brackets involving the spin variables are 
most effectively computed by considering the projection 
of the spin tensor onto the spacelike background tetrad 
vectors, i.e. S IJ = Steffi. We find 



{x\S KL }v B 
{P t ,S KL } DB 
{5 /J ,5 A ' L } DB 



oj^lu^ — 2p v UJ v 



lo^uj^ - 2p v uj u 



S ia + S^) pP (efe L a - g£g£) + 0(S 2 ) = 0{S 2 ) 



(4.25a) 



{p a uj <T ) 2 



tfzP (~L~K -L-K} +(D ( S 2) = C ( 5 2) ^ (4 25b ) 



W (etgf - etef ) (e^ - ~e%) + S IK S JL + S JL 6 



IK 



-S IL S JK -S JK S IL + 0(S 2 ) 

S IK S .JL + S JL S IK _ S IL S JK _ S JK 5 IL + ^2) 

I 



(4.25c) 



where we have used ui 1 = p 1 , which follows directly from disappear. Defining a three-dimensional spin vector by 
Eq. (|4.7p . Again the terms proportional to uj^uj 11 — 2p v u v 

(4.26) 



S I = -IJK S JK 



10 



one can immediately rewrite Eqs. (|4.25|) as 



{x\S j }ub 

{Pi,S J }uB 

{S 1 , S j }ub 



0(S 2 ), 
0(S 2 ), 
£ijkS k 



0(S 2 



(4.27a) 
(4.27b) 
(4.27c) 



Equations (|4.27[) imply that the phase-space variables 
{x\ Pj, S K } provided by the generalized NW SSC are 
canonical at linear order in the particle's spin. 



Hamiltonian in the generalized Newton- Wigner 
SSC 



In this section, we provide an explicit expression for 
the Hamiltonian (|3.27cp in the NW SSC, at linear order 
in the particle's spin. As explained in Sec. IIV Al this is 
simply obtained by inserting the NW SSC directly into 
the unconstrained Hamiltonian. Also, we express this 
constrained Hamiltonian in terms of the variables x l , Pj , 
S K , which have been proven in Sec. IIV Bl to be canonical 
at linear order in the particle spin. 

We begin by rewriting the quantity pt appearing in 
the unconstrained Hamiltonian (|3.27c[) in terms of the 
mass to = yJ—pyP^ and the spatial components pi of the 
momentum four-vector. The result is 



Pi 



-?Pi - a^/m 2 +Y 1 p I Pj, (4.28) 



where 



V^ 1 

ti 



V 



9_ 

(J* 



9 tl 9 t3 

g tt 



(4.29a) 
(4.29b) 
(4.29c) 



The crucial usefulness of Eq. (|4.28[) resides in the fact 
that the canonical phase-space variables {x l , Pj, S K } 
have vanishing Dirac brackets with the mass at linear or- 
der in the particle spin. We have established this result 
by explicit computation. As an illustration, we provide 
the details of the computation of the Dirac bracket be- 
tween x l and the mass (the other brackets involving the 
mass are computed in a similar fashion). We start from 



{x 1 , to} 



DB 



DB 



= --^—Wig^PvP^nB , 

2m 

= --^{.t\^} db , (4-30) 
to 

the last line following from {x 1 , :c j '}db = 0(S 2 ). Using 
Eq. (|3.28[) together with the fact that the Dirac bracket 
with the Hamiltonian gives the constrained equations of 



motion yields 



{x\m} DB = P^x^P^-E^S^n, 

to 

= --(p*-pV) + -/^{i'.^Idb, 
to to 

(4.31) 



where Eq. (|3.12[) must be employed to express the four- 
velocity components v l in terms of canonical variables. 
Substituting Eq. P~7) into Eq. (|3"T2"]) . it is straightfor- 
ward to show that 



i t i 
p — p V 



LO v p u 



s tx - t 

p 



(4.32) 



Next the Dirac bracket between x % and S a @ at linear 
order in spin can be computed directly following the pro- 
cedure outlined in Sec. IIVBI the result being 

W ,S^} DB = -^^(S^ + S^^) . (4.33) 
Hence, since E^ap is antisymmetric in a <-> (3, we get 



E mP {x\S a0 } 



DB 



io v p v \ u t p l ) 



UJ^p" 



u k p 
u t p l 



(4.34) 

the second line following from the definition 2E fia p = 
V AB ( eA) a {e B )i3;^ Substituting Eqs. (14321 . (14351) 
and (|4.34[) into Eq. (|4.31|) one obtains 



{x 1 , to} 



DB 



(e 1 )X;„P° 



'/A 



,txP 



+P^ T ) P ;JS^ + S^ 



: u k p 

LOtP* 



(4.35) 



Renaming dummy indices and making use of the NW 
SSC to rewrite S = —S u>k/cut, one can see that all 
terms cancel, therefore showing that the mass commutes 
with x l under the Dirac brackets. 

Since the constrained Hamiltonian depends only on 
{x l , Pj, S } and the mass to, it follows that the mass 
may be treated as a constant when taking the Dirac 
bracket between an arbitrary function of constrained 
phase-space variables and the Hamiltonian. 

Our Hamiltonian (|3.27cj) now takes the form 



H = f3 l p i + aJm 2 + ^p i p j -E tAB S AB .(4.36) 
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Equation (|4.26[) implies S 
SSC [Eqs. gl| and gTJ)] implies 

_ S"uj 



JJK qK 



s 1 



UJT 



S K . while the NW 



(4.37) 
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where 



lo t = lo^ e£ = e T - m , 



^ e i 



(4.38a) 
(4.38b) 



The canonical momenta P t are related to the linear mo- 
menta pi by Eq. (|3.20p , which may be rewritten in terms 
of the canonical spin variables as 



Pi = pi + E lAB S 



AB 



Pi + | 2E lTJ — + E lJK e JKL S L . (4.39) 

LOT 



In principle, in order to express the Hamiltonian (|4.36|) 
in terms of the canonical momenta Pi, one must invert 
Eq. (|4.39p to obtain pi as function of canonical variables 
(recall that ui^ depends on p^). However, because our 
Hamiltonian is valid only at linear order in the test- 
particle's spin, it is sufficient to write 



Pi =Pi-[ 2E lTJ — + E iJK e jn ^S^ + 0(S') , (4.40) 
where 



V. EXPLICIT HAMILTONIAN FOR SPECIFIC 
BACKGROUND SPACETIMES 

A. Spherically symmetric spacetime in isotropic 
coordinates 

The line element for a generic spherically symmetric 
spacetime in isotropic coordinates is given by 

ds 2 = -f(p) dt 2 + h(p)(dx 2 + dy 2 + dz 2 ) , (5.1) 

where p 2 = x 2 + y 2 + z 2 . The natural tetrad associated 
with this spacetime and coordinate system is 



1 

77 
i 



6$, 



pr — SM 

' Vh 1 ' 



(5.2a) 
(5.2b) 



where the symbol 5j is equal to when p = and it is 
equal to 1 when p = I numerically 12 . With a metric and 
a convenient tetrad in hand, one may now compute the 
quantity E^ab a s follows 



E 



1 



fj,AB 



[(ex)A(e^), M + (e^) A r^e Bj 



(5.3) 



P = Pi, 



(4.41a) 
(4.41b) 



Pt = -0 l Pi-a^/m 2 + i ij PiPj, (4.41c) 

(4.41d) 
(4.41e) 



Pp. &t ^ ' 



uji = ujp e^j = P^e^ . 
We may now write the constrained Hamiltonian (|4.36[) as 



H = (f Pl + a ^m 2 + 7 « Pi Pj - F t K S K + C(S 2 ) , 

(4.42) 



where 



U)T 



UK 



(4.43) 



By substituting expression (|4.40[) for pi into Eq. (|4.42|) 
and expanding to linear order in spin, one arrives at last 
at the following Hamiltonian 



H = H NS -[ /?Ff + Ff 



OK , 



sjm 2 • -'■'/',/',, 



(4.44) 



The algebra is straightforward and the result is 



E, 



fiTI 



EpjK — —-^rS^jriK] 



(5.4a) 



(5.4b) 



where the prime symbol denotes a derivative with respect 
to p, and where nj = (x/p,y/p,z/p). The last ingredi- 
ents needed in order to obtain the explicit Hamiltonian 
are lot and lok defined in Eqs. (|4.41d[) . (|4.41c|) . A quick 
computation yields 



lot = -ym 2 +j l 3PiPj 

= -m ( 1 + y/Q) , 



LO K 



P, 



K ■ 



(5.5a) 
(5.5b) 



where Q = 1 + f j P l P 3 /m 2 and P K = P 5 3 K . By substi- 
tuting Eqs. jEU) and j531) into E q- (|4~43| . we obtain the 
following expression for the quantity F^ 



1 



1 



JJK 



m(l + y/Q) Isfih 

ti 



njPj 



K 



2h 



-UK? „ 



(5.6a) 
(5.6b) 



where i?NS is the Hamiltonian for a non-spinning particle, 
simply given by 



^8 = ^ + ^+7'^ 



(4.45) 



12 More precisely, even though the spacetime index and the in- 
ternal tetrad index I are completely different in character, both 
indices may take on the same numerical value (1,2 or 3 associated 
with x,y and z respectively). 
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Finally, by substituting Eq. (|5.6[) into the Hamilto- 
nian (|4.44p and performing simple algebra, we arrive at 



H = H] 



NS 



Qif'h - fh<) - fh' 

2M P y/Th 2 y/Q(l + y/Q) 



(L-S*), (5.7) 



where -ff^s is the Hamiltonian for a non-spinning particle, 
and where 



Q 


= 1 + 




(5.8) 


p2 


. sJkPjPk _ 

9 


771 z 


(5.9) 






L-S* 


= pe IJK n lPj 


m ^ 


(5.10) 



The quantity M in Eqs. (|5.7p and (|5.10|) is introduced in 
anticipation of specialization to the Schwarzschild met- 
ric below. Since a spherically symmetric spacctime pos- 
sesses an 50(3) symmetry (associated with rotation of 
the x,y,z coordinates among themselves) that is shared 
by the internal tetrad space, one may accompany any co- 
ordinate rotation by the corresponding tetrad rotation, 
thereby preserving the functional form of the Hamilto- 
nian (|5.?p . as well as the quantities (|5.9[) and (I5.10[) . 
Thus one may meaningfully identify the vectors Lj = 
p e IJK tijPk and Si (which really live in the tetrad in- 
ternal space) with spacctime vectors Li and Si which 
transform accordingly under rotations of the coordinates 
x,y,z. 

In the limit of flat spacctime, the Hamiltonian re- 
duces to -Hns as expected, since the Cartesian compo- 
nents of the spin are all constants of motion. For the 
Schwarzschild spacctime in isotropic coordinates, we have 



ds 2 



1 ~ M/(2p) 
l + M/(2p) 



dt 2 - 



+ (dx 2 +dy 2 +dz 2 ) 
2pJ 

(5.11) 

Substituting these explicit expressions for f(p) and h(p) 
in the Hamiltonian (15.71). one finds 



H = H 



*/' h 



NS 



p 3 VQ(i + VQ) 

M / M 

1-2 1 

2p \ Ap 



(L-S*), (5.12) 



where tp = (1 + M/2p)~ 



here generally accepted notation conventions. The natu- 
ral tetrad associated with this spacetime and coordinate 
system is 





l 

— <5 M 

VT * ' 


(5.14a) 


-IX 

e l = 




(5.14b) 


e 2 = 


r 


(5.14c) 


4 = 


1 

rsinfl 


(5.14d) 



The metric (|5 . 1 3[) and the tetrad (|5.14|) then lead to the 
following result 



EtAB 
E r AB 

Eqab 



E, 



4>AB 



L XT rl 

2~m [A B] ' 

o, 

-foS[A S l] + cos6Sf A 5 3 B] 



Next the computation of lot and lok yields 



LOT 

w 3 



m ( 1 + v/Q 



1 

= -Pa 



Pr 



rsmf 



(5.15a) 
(5.15b) 
(5.15c) 

(5.15d) 



(5.16a) 
(5.16b) 

(5.16c) 
(5.16d) 



Equations (|5. 15(1 and (|5. 16[) then allow us to obtain F^. 
The result is 



Fl = cos( 



F 2 
ii 



(5.17a) 



Pi 



2rsin9y/JhJ \1 + VQ 

16 V. 



f: 



Vh 11 

-L\ ( _M * + ±* 



(5.17b) 



(5.17c) 



B. Spherically symmetric spacetime in spherical 
coordinates 



In this case, the metric takes the form 

ds 2 = -f(r)dt 2 + h(r)dr 2 + r 2 d9 2 + r 2 sin 2 Odcf 2 . (5.13) 

Note that the functions / and h appearing above are not 
the same as in the isotropic case. However we follow 



where again Pj = Pi/m. The Hamiltonian then follows 
immediately 



H = H 



NS 




COS0 p a 



(—^P4>S2 


+ PeS 3 ^j 


P^S 2 


PeS s \ 




r 2 Vh J ' 
(5.18) 
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The spin terms in the first line of the Hamiltonian (|5.18p 
are the spherical coordinate equivalent of the L ■ S* terms 
of the isotropic Hamiltonian (|5.7p . The spin terms in 
the second line of Eq. (|5. do not vanish in the flat 
space limit / = h = 1, and therefore represent coordinate 
effects related to the fact that the components of the spin 
in spherical coordinates and its associated tetrad must 
evolve, even in the absence of spin-orbit coupling. Such 
spin terms in the Hamiltonian represent therefore a type 
of gauge terms. 

Notice however that one could in principle eliminate 
these gauge terms in the Hamiltonian by picking a 
"Cartesian" tetrad, even though the coordinate system 
chosen is the spherical one. For example one could pick 
the "isotropic" tetrad (|5.2j) . taking care of transforming 
the components of &a from isotropic to spherical coordi- 
nates. In that case the spin degrees of freedom Sk, which 
live in the internal tetrad space, behave as the compo- 
nents of the spin in Cartesian coordinates, and in that 
case the flat space limit of the Hamiltonian should be free 
of gauge terms and should reduce to the non-spinning 
Hamiltonian. 

For the Schwarzschild spacctimc, / = 1/h = l — 2M/r, 
and we obtain 



H 



H 



M 



NS 



1 - 2M/r 



where 
Q = 1 



VQ 
PaS 2 



2M 
r 



r 2 sin 2 i 



,P^S 2 + P e S 3 



2M\" 1/2 - 



PeSz 



9 • 2 

r z sin 



P 2 



(5.19) 



(5.20) 



C. Kerr spacetime in Boyer-Lindquist coordinates 

Not surprisingly the computation of the Hamiltonian 
is much more involved in Kerr spacetime, whose line el- 
ement, in Boyer-Lindquist coordinates, is given by 

,2 a 



ds 2 = 



2Mr\ ,5, 4aMrsin „ 
-1 + dt 2 dtt 



A sin 2 9 



b 2 + $ dr 2 



E 
Ed6> 2 



Our choice for the reference tetrad is given by the 
"spheroidal" tetrad 



el 




2aMr sin ( 



'sm9 



(5.24a) 

(5.24b) 
(5.24c) 
(5.24d) 



which reduces to the "spherical" tetrad (|5.14p for a = 0. 
This tetrad then leads to the following components for 
the quantities E^ab 



EtTl = 



tTI — 



EtT3 
Et 12 

Et 13 
Et 23 



Mw 2 p 2 
2VAS 2 ' 
a 2 y/~AM r cos 8 sin 9 

Vas 2 

o, 

0, 

ay/AMp 2 sin 9 
27IS2 ' 
aMrw 2 cos 9 

Ie 2 ' 



E r Tl 
E r T1 



0, 
0. 



E T T3 — — 



aM (2r 2 E + n7 2 p 2 )sin 



2VAAS 



E r 12 

E r 13 
E r 23 



a cos tf sm ( 
2VAS 

0, 
0. 



(5.25a) 

(5.25b) 

(5.25c) 
(5.25d) 

(5.25e) 
(5.25f) 



(5.26a) 
(5.26b) 

(5.26c) 

(5.26d) 

(5.26e) 
(5.26f) 



(5.21) 



where 



E 
A 



A 



2 2 

a cos i 



2Mr . 



a 2 A sin 2 9. 



(5.22a) 
(5.22b) 
(5.22c) 
(5.22d) 



For sake of shortening some further formulas, we also 
introduce the quantity 



a 2 cos 2 9 . 



(5.23) 



Eoti 

EeT2 



E, 



T3 



Ee 12 

Eg 13 
E023 



0, 

0, 

fl 3 i/AAfr cos 9 sin 2 9 



AE 



Ar 



2E 
0, 
0. 



(5.27a) 
(5.27b) 

(5.27c) 

(5.27d) 

(5.27e) 
(5.27f) 
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E 



F 



aM sin 2 9 



Tl 



T2 



2\/AE 2 
a 3 VAMr cos 9 sin 3 



2r 2 E + tuV) , (5.28a) 



Vae 2 



Effj 12 : 
-^013 = 



: 0, 

o, 



(5.28b) 

(5.28c) 
(5.28d) 



A sin* 
2VAS 2 

(2Mrw 4 + AE 2 ) cos 



23 



2VAE 2 

while <Z>t and a>j<- are easily found to be 



E 2 -a 2 Mp 2 sin 2 0),(5.28e) 
(5.28f) 



m ( 1 + 
A 



LOT ~ 
LUl = P r 



^2 = Pey^ 



^3 = Pd 



1 

E' 



sin 6>V A 



(5.29a) 
(5.29b) 
(5.29c) 
(5.29d) 



where 



A - 2 1 A2 E 

=1 + -P 2 + -p e 2 + — - 

^ Asm i 



■ P 2 



(5.30) 



with p = Pijra. The coefficients F^ are finally given 
by 



Ft = 2aMrcost 



ay/A ~ w 2 

A(l + VQW 2 ~VAE^ 



F 



aM (2r 2 S + tu 2 p 2 )sin( 



VAA (1 + VQ) E 3 / 2 
2a 3 Mrcos#sin 2 6 l /~A 



Pa 



Fl = cos 9 



A(l + y/Q) V S3 
2Mrw 4 + AS 2 



Pe 



(5.31a) 
(5.31b) 

(5.31c) 



\/AE 2 



2a 3 Mr sin 2 9 I A 



A(i + VQ) V e 3 



Ai> 2 



— - 



Ft = -- 



A (1 + VQ) sin 
aM (2r 2 E + p 2 ro 2 )sin6> 

A (1 + VQ) £ 3/2 

2a 3 Mr A cos 9 sin 2 



P 



~ a 1 



P 



-sin 



P 



a(i + vqw 2 

~aM(2r 2 Z + zv 2 p 2 



(5.32a) 
(5.32b) 

(5.32c) 



2 P * 



a (i + VQ) £ 3/ 

/A ^rE 2 -a 2 Mp 2 sin 2 6i 
E 2 



(5.32d) 



P^ 

f: 



M 



VA(i + VQ)s 5/2 

P e + 2a 2 r A sin cos 9 P r ) , (5.33a) 



a" cos V sin ( 



AE 



Ar 



(5.33b) 
(5.33c) 



aM sin 2 9 



Va(i + VQ) s 5/2 

+ (2r 2 E + p 2 w 2 )P e 



2a r A cos 9 sin 6P r 



(5.33d) 



Inserting these results into the Hamiltonian (|4.44[) . a long 
but straightforward computation yields 



H — Pns + HiS 1 , 



(5.34) 



(5.31d) where 



15 



H, 



A cos (9 



A 2 VSVQ(1 + VQ) sin 2 6 
aMA(2r 2 T, + w 2 p 2 )sin9 
A 3 / 2 £ 2 VQ(1 + VQ) 



(1 + VQ)(AS 2 + 2Mrw 4 ) + 2a 2 Mrw 2 x /Qsm 2 9 
2a 3 Mr A cos 9 sin 2 9 



Pa 



P r Pe 



A3/ 2 £V0(i + VQ) 



1 



Ho = 



A(l + V^X^ 2 ~ a 2 Mp 2 sin 2 6>) - MjQ(p 2 vo± - 4a 2 Mr 3 sin 2 (9) 



A 2 VX\/Q(1 + VQ) sm9 
aA/VA(2r 2 S + w 2 p 2 ) sin 9 



Pi, 



2S p2 | A r>2 

Asin 2 ^ Ir'" 
2a 3 AfrA 3 / 2 cos6»sin 2 6»" 

A3/ 2 s 2 v^(i + v^T) . 



A3/ 2 S7Q(1 + v^) 
a 2 A cos 9 sin 6> 



2£ ~, 1 
Asm 9 T 2j 



(AE)3/ 2 VQ(1 + v/Q) 
aMVA 



A + VQAS P r - 



-AA + w 2 SVQ(^A - Af(r 2 - a 2 )) 
(AS)3/ 2 VQ(1 + VQ) 



Po 



A 2 EVQ(1 + VQ) 



2a 2 rA cos 9 sin 9P r + {2r 2 T, + w 2 p 2 )P e 



Pi, 



(5.35) 



(5.36) 



(5.37) 



Setting a = in this result and noting that for a — one coordinates [2l| 



has A = r 4 , S = r 2 and A = r(r — 2M), it is easy to 
check that this Hamiltonian reduces to the Schwarzschild 
result (|5.19[) in the non-spinning case. 



'-a 2 -f3i 
-Pi Hi J ' 



(6.1) 



VI. COMPARING THE HAMILTONIAN IN THE 
GENERALIZED NEWTON- WIGNER SSC WITH 
THE ADM CANONICAL HAMILTONIAN OF PN 
THEORY 



In this section we specialize our results to the case of 
the Kerr spacetime, but this time using ADM-TT coor- 
dinates. By expanding our Hamiltonian (|4.42[) follow- 
ing the prescription of PN theory, we verify explicitly 
that we recover the known test-particle limit results of 
the Arnowitt-Deser-Misner (ADM) canonical Hamilto- 
nian computed within PN theory alone. The latter is 
currently known through 2.5PN order for the terms lin- 
ear in the spin [l7j] , and thr oug h 3PN order for the terms 
quadratic in the spin [ij, OS S3, SI 111 SI SI- We 
cannot reproduce the PN couplings of the test particle's 
spin with itself because the MPP equations, as we have 
already stressed, are only valid to linear order in the par- 
ticle's spin. In addition we also obtain the terms linear in 
the spins at 3.5PN order of the canonical ADM Hamilto- 
nian in the test-particle limit. Those contributions have 
never been computed before. 

In order to make the PN expansion as clear as possi- 
ble, we restore factors of c in this section. However these 
factors of c must be viewed purely as dimensionlcss PN 
book-keeping parameters, and as such we are still for- 
mally employing geometric units. 

First, let us introduce the Kerr metric in ADM-TT 



-1/a 2 -/37a 2 \ , fi9 , 
-{V/a 2 7 « - U.) ' ( 6 - 2 ) 



where j zk jkj = 5* and f3 l = j zk (3k- Defining n l = x 1 jr 
and introducing a dimensionlcss three-vector x defined 
as 



"Kerr 

M 2 



(6.3) 



where M is the mass of the Kerr black hole and S^err its 
spin, the lapse function is given by |2l| 



M 1 M 2 1 M 3 

q = c h 



rc 2 r 2 c 3 4 r 3 c 5 
1M 4 l M 3 [3(x-n) 2 - X 2 ] 
8 r 4 c 7 2 r 3 c 5 

(6 . 4) 



2 r*c 
the shift vector is given by 

f 2M 2 6M 3 21 A'/ 4 



if 



r 2 (? r 3 c 5 2 r 4 c 7 
M%{ X -n) 2 - X 2 \ 



e^XjTik + O (9) ,(6.5) 



and the spatial metric j lJ is given by 

7 « = ^-<5 ifc ^ T + O(10) , (6.6) 
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where e 
r 

as 



e y is the Levi-Civita symbol (with £123 



123 _ ]^ anc j w }j ere t} ie quantities A and /i^ 1 are defined 

3(x-n) 2 ] 



.4 



1 



2rc 2 



M 3 [ X 2 



lMY 3A/ 4 (*-n) 2 
+ 2 r 4 c 8 r 4 ^ 
iTT JW, , 

n U ~ 2 r 4 c 8 ij 



(6.7) 



7 r-5 o- 



M 4 (x-n) 2 n j; n J 



+7 



21 



2 r 4 c 8 



(6.8) 



For the reference tetrad appearing in the Hamiltonian, 
we chose 



e T 



V*, 



0(8) 



(6.9a) 
(6.9b) 



It turns out however that we only need the spatial triad 
ej through order 1/c 7 for our purposes. (This makes 
the spatial triad very simple because the spatial metric 
is diagonal at that order). 

The canonical spin S 1 appearing in the Hamilto- 
nian (|4.42p scales as the physical spin of the test par- 
ticle. To conform with standard power counting in PN 
theory, this spin variable carries a power of 1/c. There- 
fore when restoring the factors of 1/c for the purpose of 
PN bookkeeping, we make the replacement 13 



<? 7 

S 1 ^ — 



(6.10) 



Finally we define the orbital angular momentum as 

/.' n. (6.11) 



and rescalcd momentum and spin as 

p = i-p, 

m 

s* = M.s, 

m 



(6.12a) 
(6.12b) 



which are useful to abbreviate formulas below. With 
these tools it is straightfroward to expand the Hamil- 



13 This is appropriate if the particle is a black hole or a rapidly 
rotating compact star. In the black hole case, 5 = am 2 /c, with 
a ranging from to 1 [see Eq. (16 - 3 D 1 . In the rapidly spinning 
star case one has S = mv ro tR ~ mcR s ~ m 2 /c (where we have 
assumed that the rotational velocity v TO t is comparable to c and 
that the stellar radius _R is of order of the Schwarzschild radius 
R s = m/c 2 ). 



tonian (|4.42|) in powers of 1/c as 



H — mc? + H N + -it iTlPN + -~ #1.5PN + —7 #2PN 

c z cr c 

+ — #2.5PN + — #3PN + — S3.SPN + O (8) 



C J C u 
^2\ 



where 



/ P 4 3M - 2 M 2 



Hi 



5PN — -3 f 2SKerr + — S* 



(6.13) 

(6.14) 

(6.15) 
(6.16) 



H 



2PN 



ml 1 P 4 + — ttP 5- 

^ 16 8r 2r 2 4r 3 J 

+ 2Mr 3 ^ nij ~ ^ ^Kcrr ("^Korr + ) ' 

(6.17) 



2.5PN 



M 
r \ 



(6S KeTT + 5S*)-^P 2 S* 



(6.18) 



^3PN 



7M p6 27M 2 pi 




-{p x s*y(p x s KeiI y 



(6.19) 



where n,- 



mrij and S£? crr 



^Kerr^Kerr- The n0n - 



spinning terms in the Hamiltonian (|6.13p coincide with 
the corresponding terms computed in PN theory in the 
test- particle limit [32l |; the linear terms in the spins at 
1.5PN and 2.5PN order agree with the terms computed 
in the test-particle limit in PN theory [l|| [I?) ; the terms 
quadratic in the spin of the larger body coincide with 
what derived in PN theory at 2PN [33|] and 3PN or- 
der [H, H3|. We nnd that the contributions at 3.5PN 
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are given by 



9m 



-ff3.5PN = o 71 r9 A (^Kcrr ' n )(S* x ^Kcrr) 1 -P 

[S^Kexr-n) 2 -^] (9S*- 

21 M 2 

4S , Kcrr) ' i H ^"T— S'Kcrr ' £ 

2r D 

, / 7 64 , 27M - 2 105 M 2 W 



16r 3 



8H 



(S*-i). 
(6.20) 



While the terms of this expression which are cubic in the 
spins (S'xorr and S^ err S*) have already been calculated 
for generic mass-ratios in Refs. [2l|, [22[, with which we 
agree in the test-particle limit, the terms linear in the 
spins are, as far as we are aware, a new result. Of course, 
because our Hamiltonian is only valid at linear order in 
the particle's spin, this result is still incomplete as it does 
not include terms (S*) 3 and SKen(S*) 2 , which are still 
unknown. 

Finally, we stress that at leading order our general- 
ized NW SSC reduces to the so-called baryonic SSC 
of Refs. [ll], Ef|- In fact, at leading order pi w mu 1 , 
mc 2 and ~ cS^, which yields oj t ~ — 2mc 2 and 
iv % . Therefore, our generalized NW SSC becomes 



Pi 

UJi 



S> 



1 ■ -v 3 

n° 1 



(6.21) 



2 c 

in agreement with Refs. [HI, EB]. 

VII. CONCLUSIONS 

In summary: starting from the Lagrangian put forward 
in Ref. 27 1 building on the classical work of Ref. [3(| on 
the relativistic spherical top dynamics, we derived the 
unconstrained Hamiltonian for a spinning test-particle 
in a curved spacetimc, at linear order in the particle's 
spin. The equations of motion for this Hamiltonian co- 
incide with the MPP equations of motion. The latter are 
well-known to describe the motion and spin-precession 
of a test-particle, but are expressed in terms of the spin 
tensor S^ v carrying six degrees of freedom. In order to 
eliminate three of these degrees of freedom (which can be 
shown to correspond to the choice of the point internal to 



the spinning body whose worldline is followed j29j|), we 
impose the so-called NW spin supplementary condition, 
suitably generalized to curved spacetime. Using the for- 
malism of Dirac brackets [3!| we computed the Hamilto- 
nian and phase-space algebra of the constrained system. 
In particular, we showed that, in a generic curved space- 
time, the resulting phase-space algebra is canonical, i.e. 
it has the standard symplctic structure for the set of dy- 
namical variables (q,p, S), at least at linear order in the 
particle's spin. As a consequence, the equations of mo- 
tion can be derived from our constrained Hamiltonian by 
means of the usual well-known Hamilton equations. 

As an application, making specific choices of the tetrad 
field, we computed explicitly the constrained Hamil- 
tonian for a spherically symmetric spacetime, both in 
isotropic and in spherical coordinates, as well as for the 
Kerr spacetimc in Boycr-Lindquist coordinates. We no- 
tice that different choices of the tetrad field would lead 
to different Hamiltonians connected by canonical trans- 
formations. Also, we expanded our Hamiltonian in PN 
orders and showed explicitly that it reduces to the test 
particle limit of the ADM canonical Hamiltonian com- 
puted in PN theory [3, M, M, M, HI- Notably, we 
recover the known spin-orbit couplings through 2.5PN 
order and the spin-spin couplings of type Snerr S through 
3PN order, Skch being the spin of the Kerr spacetime. 
Our method allows one to compute the PN Hamiltonian, 
in the test particle limit and at linear order in the par- 
ticle's spin, at any PN order, and as an application we 
computed it at 3.5PN order. 

Another application of this work will be developed in 
a follow-up paper, where we will use our Hamiltonian 
to build a new effective-one-body Hamiltonian for spin- 
ning bodies [3l|, [H, HH, [H| • Such work will be important 
to build templates for the search of gravitational waves 
with ground and space-based detectors, as it will per- 
mit taking full advantage of the analytical and numerical 
treatment of the dynamics of spinning black-hole binaries 
throughout the inspiral, merger and ringdown phases. 
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